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Global Distribution of Lightning Activity 



Mean annual global lightning flash rate (flashes knr 2 yr 1 ) derived from a combined 8 
years from April 1995 to February 2003. (Data from the NASA OTD instrument on the 
OrbView-1 satellite and the LIS instrument on the TRMM satellite.) 






Filling CMAQ Grid Cells with Flashes 

S Get N g = # CGs from NLDN 

(every hour of each day for each CMAQ cell) 

■s Get August Z ratios: 

o from 4 year 

(Boccippio et al. 2001) 
data base. 

o assigned to 3 

o (0.5° Boccippio Cell) 

o to CMAQ cells 

o 

✓ N c = ZN g = # Cloud Flashes 


Zcmoq_ouq_smthd 


VHF 

data 


Me Caul 
Clustering 
Algorithm 


Clustereii 
VHF data 


NLDN 

data 




□ Filtration of VHF Sources 

□ Flash-Typing 

□ Transformations, Spatial 
Averaging, and Sorting of 
VHF Sources 

□ Channel Length Computation 

□ Channel Segment Creation 

> location 

> polar angle 

□ NOx Computation 

> Lab [Wanget al., 1998] 

> Theory [Cooray et al. 2009] 
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NALMA Network 


LNOM 
Analysis 
Cylinder 
(vertical extent 
is 0-20km) 

Mimics a 
CMAQ grid 
volume. 

Examine all 
flashes in 
cylinder for 
August 
2005 - 2009 . 
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DC Area Lightning Discharge- Animation 


• 2.2 sec hybrid flash 

• 50 km horiz extent 

• Initiation at 5.2 km 

• VHF Sources 2187 

• CG strike at 2 s 
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Preprocessing Steps 


1. McCaul Clustering Algorithm 

s Ingests LMA data 

/ Only consider vhf sources with chi-square of 2 or less 
/ Cluster these high quality vhf sources into “ Things ” 

✓ Send output files of clustered data to the LNOM 

2. Thing Filtration 

✓ Range Filter (remove if closest vhf source outside cylinder) 

S Number Filter (remove if it has <20 vhf sources... e.g., ) 

v What remains; legitimate flashes that are at least partially in cylinder 

3. Additional VHF Source Filtration 

/ Power Filter (remove vhf source if power < Idbw) 

/ Altitude Filter (remove vhf source if at sfc or >=20km) 

4. Flash-Type Categorization 

✓ Deemed a CG if a NLDN detection is within (100 ms, 10km) 
and lowest altitude vhf is bellow N-Region. 

5. Processing of Filtered VHF Sources 

✓ Coordinate Transformation 
s Spatial Averaging 

/ Sorting 



Coordinate Transformations 



WGSS4 

Ellipsoid 


Geodetic W? ECEF 

Standard: , X g „ /O — » (X, Y,Z) 

Heikkinen (1982): (<p_ , X_ . h g ) <— (A' J'.Z) 


vhf source 


Spatial Averaging 


Must have 1 or more 
vhf sources in cube 



(100 m ) 3 
cube 


X 


nboxmin = 1 
scale = 100 m 



Channel Length Algorithm 



o Compute all the distances between each averaging point 
(each cube has 1 averaging point, unless the cube was 
ignored) 

o Store each distance, D tJ , between ith and jth averaging 
point in matrix D. 

o At this point, note that D is traceless & symmetric 

Assign any element of D that is < 0.1 m to a very large 
number 



Channel Length Algorithm (cent.) 


o Begin with n averaging points. 


o 1 st Iteration: 

- Start at the highest altitude point. Define it to be “on the 
channel”, and all other points to be “free” (off channel). 

■ Draw line from starting point to closest free point. This 
is the first channel “section”. 


> Now there are 2 points on the channel, and n-2 free points 
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Channel Length Algorithm (cont.) 


o 2 nd Iteration: 

■ Find closest free point to 1 st channel point 

■ Find closest free point to 2 nd channel point 

■ ■ 

■ Pick min of the mins (i.e., draw line from a free point to a 
channel point that is the minimum distance). 

o Continue with more iterations until no more free points 




Sample 

LNOM 

Channel 

Construction 

(Aug 2006) 


Results shown are for: 
nboxmin = 5 
scale = 1000 meters 

Recent improvements 
to LNOM now provide 
more channel resolution 
than shown here: 

->More tortuosity 
^Greater channel length 
->More NOx 
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Sample 

Channel 

Lengths 

(Aug 2006) 
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100 


Size (stats) = 8053 
Size (shown) = 7155 
Averoge = 51,438659 
Std Dev = 54,374795 
Max = 1350,7276 
Medfon = 36,204969 
Min = 1,8396797 
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Computing Segment Altitude Distribution 



Contributes 
to NOx in 
Cylinder 





>Chop up Channel into 10-m Segments 

>Tally segments in each 100 m layer 
(only within cylinder) 







(Aug 2006) 





Convert 10 m Segments to NOx 


A. Yu- Jin Wang et al. (1998) Spark Chamber Results: 

... is slightly more complicated than: 



B. Cooray et al. (2009), but generalized for LNOM: 


. . . these last 6 depend on many variables: segment pressure, segment polar angle orientation, 

# K changes/segment, dart leader speed & current, stepped leader speed & current, continuing 
current speed & current, max path length from segment to a channel termination point, various 
production coefficients, fraction of ground flashes containing continuing currents, multiplicity. 




20 


Sample 

NOx 

Profile 

(Aug 2006) 
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Huntsville @ 200 m, so 
we estimated these 2 bins - 
by downward extrapolation 
when applying to CMAQ. 
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Return Stroke 
NOx 

Contribution 
(Aug 2006) 
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Remaining NOx Contributions (Aug 2006) 





LNOM Analysis Summary Statistics 




LNOM Lightning NOx 




o Ran LNOM on 5 Augusts: Aug 2005-2009 (27,873 flashes) 
o This gives (Ground, Cloud) Flash NOx Profiles 
o Averaged the 5 Ground NOx Profiles (gives an Aug estimate) 
o Averaged the 5 Cloud NOx Profiles (gives an Aug estimate) 
o Convert NOx profiles to NOx per flash profiles 

(just divide by # of ground or cloud flashes analyzed) 

o Fill CMAQ Grid Cells with flashes (see next slide) 
o Multiply NOx per flash profiles by # flashes in CMAQ cell 
to get NOx profile in cell; 

N-Region has to be near - 15 °C (Charging Zone): 




Ground 

Flashes 

(Aug 2006) 





I Cloud 
Flashes 

N c = Z N g 


(Aug 2006) 





Impact on CMAQ 
Ozone Predictions. 
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Figure 9. Maximum increase in surface Q 3 due to lightning NQx i ppb /. 
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Summary of 2006 Analysis 


> LNOM provides most detailed “bottom-up” LNOx estimate: 

□ flash-by-flash analysis 

□ specific channel geometry/height/length accounted for 

□ Wang et al (1998) used to parameterize Return Stroke NOx production 

□ Cooray et al (2009) used to parameterize 6 additional NOx producing processes 

> Latest LNOM NOx estimates show: 

□ Ground flashes produce ~ 14X more NOx than Cloud Flashes 

□ Ground flash NOx - 484 moles but varies! 

□ Cloud flash NOx - 35 moles 


> Impact of LNOx on August 2006 CMAQ Run are Significant 




Most Important Upgrades in 2011 


>Feb 08: Changed segment length from 1 to 10 m to improve LNOM speed. 


>Mar 09: Changed nboxmin from 5 to 1, and scale from 1 km to 100 m to 
improve channel length computation. 

>Feb 23: Height of N region computed using ave monthly sea-level temp 
(not SAD). Also updated leader current values. 

> Jun 23: Added flash-specific output files. 


>Jul 15: Improved flashtype subroutine to account for ambiguous flash 

types (V2/V3 nuances & 15 kA rule). 

>Aug 04: Expand to any month (ave sea-level temps for all months in 

2005-2010). 

>Aug 18: Improved parameterization of +CG continuing currents (75% 

have, compared to only 30% of -CGs). 


>Aug 22: Simplifying continuing current source to 100 ms duration. 




Desired Future Application 



Ground Flash Fraction 
Retrieval Algorithm 


GLM = Geostationary Lightning Mapper (on GOES-R) 


Vital link 


MSFC LNOM 
(lightning NOx) 




# GLM 

Ground 

Flashes 


# GLM 

Cloud 

Flashes 


Air Quality Models (e.g. CMAQ) 

Global Chemistry/Climate Models 

(e.g. GISS Model E, Geos Chem) 







More details on Ground Flash Fraction Retrieval 


Koshak, W. J., Optical Characteristics of OTD Flashes and the 
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2 7, 1822-1838, 2010. 

Koshak, W. J., R. J. Solakiewicz, Retrieving the Fraction of Ground 
Flashes from Satellite Lightning Imager Data Using CONUS-Based 
Optical Statistics, J. Atmos. Oceanic Technol., 28, 459-473, 2011. 

Koshak, W. J., A Mixed Exponential Distribution Model for Retrieving 
Ground Flash Fraction from Satellite Lightning Imager Data, J. Atmos. 
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